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Nonlocal entanglement is crucial for quantum information processes. While nonlocal entanglement 
has been realized for photons, it is much more difficult to demonstrate for elec tron s. One approach 
that has been proposed is to use hybrid superconducting/normal- metal devices^'^', where a Cooper 
pair splits into spin-entangled electrons in two spatially separated normal-metal leads. This process of 
nonlocal Andreev refle ction is predicted to lead to a negative nonlocal resistance and positive current- 
current correlatiorF*^. By cross-correlation measurements as well as measurements of the local and 
nonlocal resistance, we present here experimental evidence showing that by independently controlling 
the energy of electrons at the superconductor /normal- metal interfaces, nonlocal Andreev reflection, 
the signature of spin-entanglement, can be maximized. 

Multiple particle entanglement has been realized for photons^ and has been used for quantum cryptography. How- 
ever, for electrons, which are massive fermions, nonlocal entanglement has not been clearly demonstrated. Entangle- 
ment for electrons may arise in the spatial degree of freedom (orbital entanglement) or the spin degree of freedom (spin 
entanglement). Recently, orbital entanglement in a fermion ic Ha nbury Brown and Twiss two-particle interferometer 
was observed using current cross-correlation measurement^^*^, but further investigation is still required to verify 
whether the measured correlation is truely due to entangled stated. Spin entanglement has been predicted to exist 
in hybrid superconducting/normal-metal (SN) devices^ ^ . In this case, a Cooper pair splits into two spin -ent angled 
electrons in spatially separated normal-metal leads in a process called nonlocal Andreev Reflection (AR^PEI, which 
leads to a positive current-current correlation and negative nonlocal resistance. Although a cross-correlation measure- 
ment in these systems would pro vide direct evidence of entanglemenifiSlj previous experimental attempts focused on 
nonlocal resistance measurementd^^ ^^ l ^'^ l ^^ l Here we show evidence of nonlocal entanglement from cross-correlation 
measurements as well as measurements of the local and nonlocal resistance. By independently controlling the energy 
of electrons at the superconductor/normal- metal interfaces, we show that nonlocal Andreev reflection, the signature 
of spin-entanglement, can be maximized. 

Spin-entanglement at the SN interface can be understood in the context of ABP^, in which a low energy electron 
in the normal metal impinges on the SN interface and a hole is retroreflected while a Cooper pair is created in the 
superconductor. When two normal metals are coupled to a superconductor with spatial separation comparable to the 
superconducting coherence length (^s), roughly the size of a Cooper pair, it is predicted that electrons i n th e two 
normal metals can also be entangled via a nonlocal analog of AR called crossed Andreev reflection (CAR)2M1. As a 
Cooper pair splits into two entangled electrons that are then injected into the two normal-metal leads, instantaneous 
currents of the same sign are generated across the two SN interfaces, giving rise to a negative nonlocal resistance as 
well as a positive current-current correlation between the SN junctions. 

The nonlocal resistance of hybrid SN devices has been intensively investigatecP^^^'^^'i^, but the observation of CAR 
is complicated by another nonlocal process called elastic cotunneling (EC), in which electrons in the normal-metal 
leads tunnel across the superconductor with the help of Cooper pairs, resulting in a positive nonlocal resistance and 
a negative current-current correlation. Theoretical studies have found that the two nonlocal processes tend to cancel 
each other exactly in the lowest order approximation in the tunneling limili^, but experimental results and further 
theoretical investigation show that exac t canc ellation does not occur if the normal-met al lea ds are ferromagneticJ^^', if 
the NS interfaces are highly transpareni! ^^ ! ^^ !, or if there is strong Coulomb interactioiJ^^Ell. F or tun neling junctions 
the Coulomb interaction may lead to a transition from EC to CAR with increasing voltage bia^^^Ell 

The noise cross-correlation measurement is an especially powerful tool to probe the correlation between charge 
carriers in mesoscopic systemJ^^l Por example, cross-correlation measurements in the Hanbury Brown and Twiss 
type experiment have demonstrated the Fermionic nature of electrond^^EI] and have been used for probing orbit 
entanglement of electrons from independent source^^. It was also predicted that current-current correlation can be 
used to directly probe CAR and EC without the drawbacks of the resistance measurement^*^ However, despite 
considerable amount of theoretical work, no experimental observation has been reported. Below we describe both 
nonlocal resistance and noise cross-correlation measurements, and a clear signature of entanglement is observed in 
each type of measurement. 
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FIG. 1: Schematic of the cross-correlation measurement. Scanning electron micrograph with false color enhancement of a 
three-terminal superconductor/normal- metal device. The scale bar is 1 /j,m. The area colored in dark purple is Al film, and the 
area colored in gold is Cu film. The thin white stripes are due to Al film deposited on the sidewall of the lithography mask. 
The two Al/A10a;/Cu junctions are each biased with a dc current source and the ac voltage fluctuations across each junction 
are measured simultanerously. The right inset shows a schematic of the crossed Andreev reflection process. 

In Fig. 1 we show a scanning micrograph of a hybrid device, as well as the schematic of the measurement circuit (see 
the Methods section for sample fabrication and measurements). Two normal-metal/insulator/superconductor (NIS) 
junctions are formed between Al and Cu leads, with room temperature resistances in the range of 10-20 kfi. The 
junction size is about 0.3 /im by 0.45 /im, and the distance between the two junctions is about 0.28 /zm, comparable 
to the superconducting coherence length of diffusive Al. Three batches of devices were measured and consistent 
results were obtained. Here we concentrate on a single device with the most complete data. The advantage of using 
tunneling junctions instead of junctions with highly transparent interface is the following: First, the large tunneling 
resistance ensures there is little leakage current from one junction to the other, i.e., currents from different sides only 
flow to the superconducting lead (ground), so the bias across the two junctions can be varied independently. Second, 
for a highly tr anspar ent interface, EC is expected to dominate over CAR at all tempe ratures and at energies below 
the gap voltag^i^'i^, while for tunneling junctions CAR can dominate EC at high biatPEEl. Third, the Coulomb 
interaction between the electrons in the normal metals may lead to a "dynamical" Coulomb blockade effect (DCB), 
which prevents both electrons from being injected into the same normal-metal lead and thus facilitates the separation 
of the two spin-entangled electrons^. Note that this is different from the DCB in the superconducting lead which 
favors EC instead of CARi^. Fourth, tunneling junctions can be treated exactly in perturbation theory, which makes 
the comparison between theoretical predictions and experimental results straightforward. Finally, the fact that the 
junction resistance is much larger than that of the normal-metal lead ensures a sensitive measurement of current 
across the junction, as voltage is the quantity actually measured here. 

The properties of our NIS junction, as shown by the current voltage characteristics (CVCs) at different temperatures 
in Fig. 2a, can be understood by a simple semiconductor modeP^l. Below the superconducting transition temperature, 
the number of states in the superconductor available for quasiparticle tunneling decreases exponentially, leading to a 
diminishing current at low bias voltage. The ratio between the differential resistance at zero bias and the normal state 
junction resistance is (\/A:BT/27rA)e'^/'°-^"^ in the limit fc^T <C A. If there is no pin hole in the thin insulating layer, 
current flowing across the NIS junction consists of well separated events of quantum mechanical tunneling of electrons 
(no temporal correlation and inelastic scattering) . In this case, the power of the current fluctuation is proportional 
to the magnitude of the current. This so-called shot noise reflects the particle nature of the charge carrier&iSi Here 
the voltage noise power Sy is measured and it is related to the current noise power 5/ by Sy — Sj ■ (^)^, with 

whe re dV /dl is the differential resistance at bias voltage V , and F is the Fano factor, which equals 1 for full shot 
nois«Pl221. At eV 3> fcsT, the hyperbolic cotangent term approaches 1 and Si is proportional to /, the hallmark of 
shot noise, as shown in Fig. 2b for our sample. Below the superconducting transition temperature, Andreev reflection 
may lead to doubled shot noise near zero bias, as has been observed for short diffusive NS junctions with transparent 
intcrface^2ll^ for which the differential resistance at zero bias is sma ller than the normal resistance. However, for 
typical NIS junctions like ours, the probability of AR is very smalP^^^H^ thns the subgap current due to AR (the 
only charge transfer process at zero temperature) is also small. As the shot noise is proportional to the current, it 
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FIG. 2: Junction characterization, a, Current voltage characteristics of an SN junction at bath temperatures: 0.1 K (cyan), 
0.2 K (red), 0.3 K (green), and 0.4 K (blue). Inset: Bias dependence of the differential resistance measured at 0.4 K (note the 
log scale of the y axis). The solid lines are fits using the semiconductor model, with slightly elevated electron temperature, b. 
Bias dependence of the voltage noise power of an SN junction at 4.2 K (red) and 0.4 K (blue). Inset: Zoomed in regime for 
the 0.4 K data near zero bias. The solid lines are fits to the shot noise expression (Fano factor equals to 1) with background 
noise from the measurement setup also taken into account. The diverging behavior near zero bias is due to the large differential 
resistance as shown in a. 



is difficult to verify the doubled shot noise for NIS junctiond^. Although CAR can also be inferred from shot noise 
(auto-correlation) measurement^^, observing CAR by shot noise measurements is even harder since CAR is usually 
dominated by the AR process, and its amplitude decays exponentially with the distance between the two normal 
metals^" Another practical issue for auto-correlation measurements is that the extrinsic current noise from the 
measurement setup (back action) can be much larger than the intrinsic shot noise. 

Fortunately, with the nonlocal configuration shown in Fig. 1, the cross-correlation between Va and Vb is only 
determined by nonlocal processes, independent of shot noise of both junctions. In addition, it should not be affected 
by any extrinsic current noise, as there is no circuit segment shared between the two channels except a short piece 
of superconducting Al wire, along which no voltage drop is expected to occur. Moreover, the large differential 
resistance near zero bias (see Fig. 2) leads to a large nonlocal voltage-voltage correlation signal which greatly helps 
the observation of the small nonlocal current-current correlation. However, the large differential resistance here also 
limits the bandwidth used for the cross-correlation measurement due to the capacitance of the electrical wiring (see 
the Methods section), which prevents us from making a reliable measurement below about 0.2 K. 

Figure 3 shows the voltage noise power measured by cross-correlation, Sy, for our device at three different temper- 
atures: 0.4 K, 0.3 K, and 0.25 K. For SN tunneling junctions in the abscence of the Coulomb interaction in the leads, 
it was shown that when e\V\, ksT ^ A, the current-current correlation between junction A and junction B id^ 
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where Gq = 2e^//i, A^^^{A^'^) is the amplitude of CAR(EC), determined by the distance between the two SN 
junctions and the properties of the interfaces. For simple tunneling junctions, A'-'^^ — A^^ is predictecP^. Since 
the precise value of A'^^^ is not known, a normalized Sab is plotted in Fig. 3 to compare with experimental results. 
It is easy to see that by setting Va — Vb, the EC term in the bracket of Eq.(2) reduces to a bias independent 
term 2kBTA^^, and the CAR term is maximized. Since current sources instead of voltage sources were used in the 
experiment, current voltage characteristics from the semiconductor model (see Fig. 2) were used to convert Si{Va, Vb) 
to Sv{Ia, Ib) for comparison. 

At all three temperatures, when the bias of the two junctions is of the same polarity, there is a clearly positive 
correlation, a signature of entanglement, consistent with the theoretical prediction. When the bias of the two junctions 
is of opposite polarity, Eq. (2) predicts a negative Sy- In the experiment, Sy changes from slightly positive to slightly 
negative as temperature decreases, again consistent with predictions, but also suggesting a slowly increasing amplitude 
of the EC process, which means the assumption pp^^ — A^'^ at finite temperature is oversimplified. At Va, Vb 0, 
while Sab is predicted to be 2kBT{A'^'^^ — A^'^'), the experimental data show a positive correlation at 0.4 K that 
evolves to a sharp negative dip at 0.25 K. Such a sharp negative correlation again indicates that the EC process 
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FIG. 3: Bias current dependence of voltage noise power by the cross-correlation measurement. The left column shows the 
measured voltage noise power at three different bath temperatures: 0.4 K (a), 0.3 K (c), 0.25 K (e). The right column shows 
the respective theoretical prediction at each temperature, with the assumption A^^^ — A^^' in Eq.(2). 



strongly dominates over the CAR process at lower temperature, which is not consistent with theoretical predications. 

To understand this discrepancy. Coulomb interaction in the leads need to be considered ^ * ^^ * ^^ !. We note that 
the charging effect of solitary junctions has a strong temperature and bias dependencd^. As the Coulomb energy 
associated with charging the superconducting lead and charging the normal-metal leads changes, A^^^ and A^^ also 
change and this may lead to a nonmonotonic temperature and bias dependence of Sab- Note that here the junction 
capacitance is much smaller than that of the planar NISIN devices studied befor^^^l^ Since the charging energy is 
e^/2C, where C is the capacitance of the tunneling junction, a much stronger DCB effect is expected. The effect 
of DCB on noise correlation was recently considered with a renormalization group approachPSl^ where the presence 
of Coulomb interaction in the normal-metal leads is treated as a reduction of interface transparency, resulting in a 
change of the sign and magnitude of the correlation. Although this approach is more general than the analytical 
result in Eq. (2), it does not seem to be applicable as the Coulomb interaction in the superconducting lead was not 
taken into account in the theory. Further experiments with a modified local impedance near the normal-metal leads 
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FIG. 4: Bias current dependence of local and nonlocal differential resistance. The left column shows the local differential 
resistance at bath temperatures 0.4 K (a) and 0.3 K (c), obtained by adding a small ac current modulation to I a and 
measuring the ac voltage modulation of Va- The nonlocal resistance is obtained at the same time by measuring the voltage 
modulation of Vb, as shown in the right column at 0.4 K (b) and 0.3 K (d). The black dashed lines indicates a strong reduction 
of the local differential resistance, approximately where the dc bias voltage is the same for the two junctions (the normal state 
resistances of the two junctions are 15 and 19 kfi respectively). 



that enhances the DCB effect may facilitate comparison with theorjEHIl note also that in most theoretical work 
only pair tunneling (AR, CAR, and higher order processes) is considered in the zero temperature limit, while in 
the experiment, quasiparticle tunneling dominates at finite temperatures since a reasonable differential resistance is 
required for measurement. 

The local and nonlocal differential resistance corroborate the observation from the noise cross-correlation measure- 
ments: They also demonstrate a clear signature of entanglement, as shown in Fig. 4. As noted before, since the 
tunneling junction resistance is large, there is little current redistribution (flow from one normal-metal lead to the 
other via the superconducting lead), which enables us to do a 2D scan of both local and nonlocal differential resistance, 
in contrast to previous ID nonlocal resistance measurementJSnHIIIIll. Surprisingly, even for the local resistance, there 
is clearly a reduction of resistance when V^'^ « V^'^, as indicated by the dashed line in Fig. 4. This strong reduction 
is not expected as the differential resistance is determined by quasiparticle tunne ling, and even if there is a small 
probability of pair tunneling, direct AR is expected to dominate CAR and EC! ^° * ^^ l However, this might suggest that 
DCB indeed enhances the splitting of a Cooper pair into separate normal-metal leads, especially when the electrons 
passing the two SN interfaces are of the same energy, in agreement with predictionJ^!. 

Figure 4b shows that at 0.4 K, the nonlocal resistance is negative in the subgap regime for both junctions, a clear 
signature of spin entanglement. As the temperature decreases to 0.3 K, a peak evolves near zero bias (at the peak, 
the nonlocal resistance even becomes positive at 0.25 K, data not included here). Remarkably, if we take the ID cut 
a,t Ib — 0, the bias dependence of the nonlocal resistance is close to that previously reported where only one junction 
was biasecP^"^^. This confirms that the nonlocal resistance measured in the current setup has the same origin; though 
a more completed picture can be obtained here. 

In summary, we have shown experimental evidence of spin entanglement revealed by noise cross-correlation mea- 
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surements between two spatially separated normal-metal leads. We find a clear positive correlation when the bias 
of the two junctions is of the same polarity, indicating the nonlocal Andreev reflection and spin entanglement of 
electrons is maximized, consistent with theoretical predictions. At lower temperature and near zero bias, we find 
a strongly negative correlation, possibly due to a dynamical Coulomb blockade effect in the superconducting lead. 
We also find that both local and nonlocal differential resistance measurements can be used to demonstrate that the 
entanglement is maximized with the same bias voltage across the two SN interfaces. These findings may lead to better 
understanding and control of quantum entanglement devices made of hybrid SN structures. Further experiments with 
spin selective normal- metal leads and Bell inequality type measurements^^ will provide a more conclusive check of 
spin entanglement. 

Methods 

The devices were fabricated using standard two angle electron-beam lithography and e-gun evaporation. A poly- 
methyl methacrylate/polydimethyl glutarimide (PMMA/PMGI) bilayer was spin-coated on Si substrates with 300 
nm Si02 insulating layer for patterning the devices. For the particular device reported here, a 27 nm 99.999% pure 
Al film was deposited first in an e-gun evaporator with a base pressure of 2.2xl0~^ Torr at a 40 degree angle and at 
a rate of 0.1 nm/sec. After deposition, 0.2 Torr pure oxygen gas was allowed into the chamber for about 5 minutes 
to create a thin layer of oxide. Then a 50 nm of Cu film was deposited at a normal angle. Devices from three batches 
were measured in an Oxford KelvinOx 100 or an Oxford KelvinOx 300 dilution refrigerator. Results obtained on 
these devices were consistent with each other. 

The current sources shown in Fig. 1 were realized by putting a large ballast resistor (1 Gf2) in series with a filtered 
voltage source. The ac voltage signals from the junctions were amplified using home-made battery-powered low noise 
amplifiers. The voltage signals after amplification were sent to data acquisition cards and were analyzed by a computer 
after digitization. To reduce noise coupled to devices at low temperature, 7r-filters on top of the cryostat were used. 
The measurement signals from the cryostat were amplified inside a /i- metal enclosure close to the cryostat. The 
bandwidth of measurement was limited by the roll-off effect due to the capacitance of the 7r-filter (about a few nF). 

The auto-correlation data shown in Fig. 2b at 4.2 K were measured in the frequency range from 500 to 1000 Hz, 
while at 0.4 K different frequency ranges, from 200 to 600 Hz and from 120 to 130 Hz, were used and similar results 
were obtained. The voltage noise power Sv is fitted to, 

S. = 2./cotl,(^,(f)" + 2,/,*)'H-V,?, (3) 

where /„ is the extrinsic current noise, and Vn is the extrinsic correlated voltage noise. A good fit with data was 
found by using /„ =30 iA/^/Hz, and Vn =0.7 vN /^J Hz. These values are consistent with standard JFET amplifiers. 
Hence, in the sub-gap regime, the extrinsic current noise dominates the intrinsic shot noise. 

The cross-correlation data shown in Fig. 3 were obtained in the frequency range from 2 to 6 Hz to avoid the roll-off 
effect and the side lobe of the 60 Hz line frequency peak. The efficiency of the cross-correlation measurement at this 
frequency range was compromised by the 1/f noise of amplifiers. Nevertheless, due to the large differential resistance 
the bias dependence of Sab is clearly visible. 
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